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ABSTRACT 

The stellar mass-luminosity relation (MLR) is one of the most famous empirical "laws", 
discovered in the beginning of the 20th century. MLR is still used to estimate stellar masses 
for nearby stars, particularly for those that are not binary systems, hence the mass cannot be 
derived directly from the observations. It's well known that the MLR has a statistical dispersion 
which cannot be explained exclusively due to the observational errors in luminosity (or mass) . It 
is an intrinsic dispersion caused by the differences in age and chemical composition from star to 
star. In this work we discuss the impact of age and metallicity on the MLR. Using the recent data 
on mass, luminosity, metallicity, and age for 26 FGK stars (all members of binary systems, with 
observational mass-errors < 3%), including the Sun, we derive the MLR taking into account, 
separately, mass-luminosity, mass-luminosity-metallicity, and mass-luminosity-metallicity-age. 
Our results show that the inclusion of age and metallicity in the MLR, for FGK stars, improves 
the individual mass estimation by 5% to 15%. 



1. Introduction 

Mass is a fundamental stellar parameter with a 
crucial impact on the evolution and the internal 
structure of stars, and in other astrophysical sit- 
uations (e.g. the initial mass function, the mass- 
luminosity relation, the study of extra-solar plan- 
ets, the definition of the star-brown dwarf limit, 
etc) necessary to study stellar populations. Stellar 
masses can be accurately determined only for vi- 
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sual binaries with known orbital elements and dis- 
tances, for detached double-lined eclipsing bina- 
ries, and for resolved spectroscopic binaries. The 
latter can provide even accuracies of 1% to 3% 
( Torres et al.|2010 |. Unfortunately, the number of 
stars for which accurate stellar masses are avail- 



able is less than 200 


Torres et al. 


2010 


Hillen- 


brand and White|2004 


. Therefore, the estimation 



of mass for single stars makes use of indirect meth- 
ods: either observational — as the use of surface 
gravity obtained from the detailed spectroscopic 
analysis (Santos et al. 2004) — or by means of 



comparison between observations and theoretical 



stellar models (Holmberg et al. 2009 1 or asteroseis- 
mology ( Metcalfe et al.||2012 ). However, on those 
estimations the mass-luminosity relation (MLR) 
is particularly useful as, in principle, it is purely 
observational (hence, model-independent) and re- 



quires only the knowledge of luminosity (easily de- 
termined for nearby stars). This fundamental re- 
lation was definitively established by [Eddington 
(1924) triggered by previous works (for the his- 
torical MLR overview, please see Lecchini 20071. 



During the last 20 years, thanks to the observa- 
tional improvements (mainly in parallax and in- 
frared photometry) several works have been pub- 
lished on the MLR subject. These works have in 
common the derivation of formulae accounting for 
the mass-luminosity dependence for a large range 
of masses: from the brown dwarf mass limit until 
O and A spectral type stars. We can find poly- 
nomial fits between the mass and luminosity loga- 



rithm (or the visual magnitude) to be linear ( Ces- 



ter et al. 19831, parabolic (Henry and McCarthy 



1993} , cubic or higher (|Reid et al. 2002) 
overview, see 



For an 



Xia and Fu (2010). Other than the 



derived monotonic relations, all the authors point 
out that the observational MLR shows a dispersion 
that cannot be explained by the observational er- 
rors neither in mass nor luminosity. This is well 



summarized by Torres et al. (2010): The signifi 



cant effects of both stellar evolution and abundance 
differences are well seen in the close-up of the de- 
ceptively tight mass-luminosity relation of Fig. 5 
that we show in Fig. 6. Making the error bars 
visible highlights the fact that the scatter is highly 
significant and not due to observational uncertain- 
ties. In spite of all the works carried out during 
the last years on the analysis of age and chem- 
ical composition contribution to the MLR (e.g. 



Hurley et al. 20001, a MLR involving simultane- 



ously mass, age, chemical composition, and lumi- 
nosity was never established. As far as we know, 



only Bonfils et al. (20051 quoted empirical formu- 



lae where the visual magnitude is dependent both 
on mass and on metallicity ([Fe/H]) but just for M 
stars. The age was not considered. We point out 
also a kind of mass-luminosity relation derived by 



Torres et al.| (|2010j) where they performed a fit to 
mass expressed by a polynomial in effective tem- 
perature, gravity, and metallicity. However, in this 
case, a ge and luminosity are not explicitly used. 
The main goal of our work is to discuss the metal- 
licity and age contribution to the MLR. For that, 
we use the data for FGK stars from lTorres et al.l 



(2010) for which individual values of stellar mass, 



luminosity, metallicity, and age are available: this 
gives 13 binary systems. The mass estimations of 



these binary components have an accuracy bet- 
ter than 3%. This work is organized as follows: 
i) in Section [2j we derive the empirical relation 
for the mass as function of luminosity, metallicity, 
and age; ii) in Section [3] we discuss the results and 
draw the conclusions. 



Mass-luminosity-metallicity-age 
tion: fits and results 



rela- 



Both theory and observations indicate that the 
Sun was 20% to 30% fainter in ZA MS than today 
(Sackmann and Boothroyd 2003). However, the 



solar mass remains basically the same during the 
solar evolution. Therefore: same mass, different 
luminosity. The impact of chemical composition in 
the MLR is less clear. But, it can be shown, using 
simple stellar homological relations, that metal- 
licity can have an impact of O.25M for solar-like 
stars (see Section]!]). Consequently, we aim to dis- 
cuss the impact of metallicity and age on the MLR. 
It is well known that age has a strong impact on 
the luminosity of a star. In order to accomplish 
our work, we have selected 13 binary systems with 
known mass, luminosity, metallicity, and age for 
each binary member, taken from the recent com- 



pilation of Torres et al. (2010), where components 



are all FGK main sequence stars: the mass ranges 
from 0.8 to 1A5Mq; the metallicity from -0.6 to 
0.4 dex; and the age from 0.6 to 7.0 Gyr. The 
choice to restrict our analysis to FGK main se- 
quences stars lies on the age determinations. Stel- 
lar age is currently determined by means of the- 
oretical models and/or isochrones. These stellar 
models are based on the input physics to describe 
the evolution of the Sun. This approach seems to 
be suitable to reproduce FGK stars, because their 
global parameters are well recovered by means of 



theoretical models (see Lebreton et al. 1999 2008 1 . 
So, we expect that the ages of our sample are reli- 
able. As the errors on age are not given by | Torres] 
etaD(2010| we assume a typical error of 20%. Ta- 
ble^ shows the observations of the selected stars. 
The Sun is also shown. 

Using the observational data on Table[T]we per- 
form multi-dimensional fits, using inverse prob- 
lem techniques associated to least square methods. 
Such guarantees a minimum deviation between the 
considered model considered and the observational 
data points ( |Menke|[l989| . The matrixial form of 



2 



the equation we have to solve is: 

\ogm obs = Gd, (1) 

where m obs is the matrix with observed masses, 
G is the matrix that contains the information re- 
lated to luminosity, metallicity, and age, and d is 
the matrix of parameters that we obtain using in- 
verse problem techniques. 

The empirical model which we will consider is lin- 
ear in the parameters and goes up to the third 
degree of the input parameters, as presented in 
equation [2j 

N 3 

logm = ^2^2a itj xl (2) 

i=l 3=1 

Where a, are the constant parameters, x describe 
the observational values, the index j translate the 
exponent of x. The index i describe the cases 
where Log L/L Q is considered (i = 1), the case 
where metallicity is considered (i = 2) and, finally, 
when Stellar Age / Age & is used {i = 3). Depend- 
ing on how many parameters are used on the fit: 
N can be equal to 1 when only luminosities are 
being used; can be equal to 2 when both luminos- 
ity and metallicity are being used; and, finally , 
equal to 3 when all elements are being used. 

Because our matrix system is overdetermined 
we have to rewrite equation [l] to determine the 
parameter in d. We obtain the final equation con- 
sidering the least square approach shown in the 
equation [3j 

d est = {G T G)' 1 G T logm obs (3) 

To take into account the errors in the measure- 
ments we use Monte Carlo techniques, in order to 
estimate how these influence the final results on 
the fit and on the calculated mass. The equation 
[3] is computed 10000 times using gaussian random 
simulations of data points. Each simulated data 
point has the standard deviation associated to its 
observational error (see Table [lj. Note that when 
N = 1 the Monte Carlo estimated error bars are 
very small, because the errors in the observed lu- 
minosity are also small (compared to the absolute 
value of brightness) and they do not significantly 
change the final adjustment. Naturally the indi- 
vidual errors are larger for N — 2 and N = 3. 



Equations |1J [5] and [6] show the relation we ob- 
tain for N = 1, with R 2 = 0.839, N = 2, with 
R 2 = 0.898, and N = 3, with R 2 = 0.955, respec- 
tively: 

\ogm obs = 0.208(±0.011) J Lo ffj L/ J L Q 

+0.063(±0.007)(Lo 5 L/L Q ) 2 
-O.13O(±O.O48)(£o 5 L/L ) 3 (4) 



logm obs = 0.238(±0.014) J Log J L/ J L o 

+O.13O(±O.O41)(Lo 5 i/L ) 2 
-0.280(±0.081)(Lo 5 i/L Q ) 3 
+0.044(±0.027)[Fe/ff] 
-0.167(±0.075)[Fe/ff] 2 
-0.116(±0.164)[Fe/#] 3 (5) 



logm obs = O.219(±O.O23)Lo 5 L/L 

+O.O63(±O.O6O)(Lo ffj L/L ) 2 

-O.119(±O.112)(Lo 3j L/L ) 3 

+0.079(±0.031)[ J Fe/i/] 

-0.122(±0.119)[ J Fe/i/] 2 

-0.145(±0.234)[ J Fe/i7] 3 

+0.1U(±0.062)StellarAge/AgeQ 

-0.224(±0.l04)(StellarAge/Age Q ) 2 

-0.076(±Q.045)(Stellar Age/Age^ 

On Figures [T] [2j and [3] we plot the calculated 
mass using the parameters obtained by the fit as a 
function of observed mass. To help with the visual 
interpretation of the quality of the model we draw 
the one-to-one line on each plot. 

3. Discussion and conclusion 

In the previous section we derive three MLR: i) 
the classical MLR where luminosity and mass ap- 
pear explicitly; ii) a MLR including the age con- 
tribution; and hi) the MLR including both metal- 
licity and age contributions. From Figures [l] [2] 
and [3] we can see that our fits reproduce rather 
well the observed mass — solar mass is exactly re- 
produced in Figures [T] and [2j and for Figure [3] we 
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find 0.98-Mq) — sowing a natural increasing devi- 
ation for higher masses and for more evolved stars. 
However, the quality of the fit improves from Fig- 
ure [T] to Figure |3j which can be seen by the in- 
crease of the correlation coefficient from ~ 0.84 to 
~ 0.96 R On the other hand, our results shows 
(Figure |4| that if the fit uses only the luminosity, 
10 stars among the total of 26 have a predicted 
mass outside the range m obs ± O.O55Af (typical 
uncertainties for MLR kind relations |Henry and| 
McCarthy|1993 e.g. ); when using both the lumi- 



nosity and the metallicity, only 6 stars miss that 
range; and when fitting using luminosity, metal- 
licity, and age the number of bad predictions is 
reduced to 3 stars. These results show that the 
inclusion of both the metallicity and the age on 
the MLR improves the mass predictions. 

We estimate that, on average, the non-inclusion 
of metallicity and age on the mass predictions 
through the MLR, can over/under estimate the 
stellar mass by O.O5M0 (with a maximum value 
of about O.15M ) 

4. Annex: Chemical composition impact 
on the MLR. Homological approach. 

According the homology approach applied to 
the stellar structure equations and assuming e = 
eop x T u and K = n p n T~ s for the energy produc- 



tion rate and the opacity, respectively (c/. Cox 
and Giuli||T968"l ), we can write: 



4d(lni?)-nd(lnp)+(4+s)d(lnT)-d(lnL) = d(lnM) 

(7) 

On the other hand re-writing the above equation 
in relation to Sun, we can derive: 



L Q (r) 
where 



£00 



( _«o\ P ( Mo V 
\k qJ VMgo/ 



M 



(8) 



3A 



3\ + v 

h V — S - 

3A 



3n 
- v 



7 = v 



3A + v — s + 'in 
(u-s- 4) (3A + v) 
3A + v — s + 3n 



Note that the correlation coefficient is a non-linear param- 
eter ranging from -1 to 1. That is, the difference between 
0.8 and 0.9 is much higher than the difference between 0.1 
and 0.2, for instance. 



5 = A 



1 - 



(3A + v) (A + v - s + n - 2) 



3A + v — s + 3n 

We want this equation written explicitly as a func- 
tion of the chemical composition (X, Y and Z). For 
that we consider that the main source of opacity 
comes from the bound-free and free-free contribu- 
tions: K — + Kff, where 

K fc/ cx4x 10 25 Z (1 + X) pT~ 

K// « 4 x 10 22 (X + Y) (1 + X) P T- 

(10) 

On the other hand assuming the pp-chain 



3 5 cm 2 g 1 



~ 3 - 5 cm 2 g 



(9) 
-l 



e{pp) 



2.4 x 10 4 pX 2 

2 

rp 3 
1 9 



and Y ~ 2Z + Y p (where Y p ~ 0.25 is the primor- 



dial helium value, Tsivilev 



( |2009[ )) we can finally 



derive a MLR depending on metal abundance Z: 



/ M 



5. -id 



l-3Z-Y p 
1 — 3Zq — Y p 



26.00 



(2-3Z- Y p )(l-Z + 10 3 Z) 



(2 



3Z Q 
2 - 



-Y P )(1 
3Z G ~ 



Z, 



by 



2-3Z- %Y V 



f 10 3 Z, 

-7.77 

X 



0J 



L(r) 
L Q {r) 



(11) 



Therefore, using this simple approximation, we 
can conclude that for the solar luminosity a change 
of the metal abundance (Z) from 0.004 and 0.03 
(typical for the metallicity ranges of this work) 
can have an impact of about 0.25M Q . This value 
is luminosity independent for FGK main sequence 
stars. 
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Fig. 1. — Comparison between observed mass and 
mass obtained by the model when luminosity is 
considered. The one-to-one solid line is plotted to 
guide the eye. 
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Fig. 3. — Comparison between observed mass 
and mass obtained by the model when luminosity, 
metallicity, and age are considered. The one-to- 
one solid line is plotted to guide the eye. 
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Fig. 2. — Comparison between observed mass and 
mass obtained by the model when luminosity and 
metallicity are considered. The one-to-one solid 
line is plotted to guide the eye. 



Fig. 4. — The predicted mass minus the observed 
mass value as function of the observed mass: mass- 
luminosity relation (diamond), mass-luminosity- 
metallicity relation (star) and mass-luminosity- 
metallicity-age relation (cross). The horizontal 
line defines the ±0.055M Q range. 
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Table 1: Observational data from Torres et al. (2010) 



Name 


Log L/Lq 


error 


[Fe/Hl 


error 


Stellar Aqe / Aqep) 


error 


Mass 


error 


V570_Per A 


0.6580 


0.0230 


0.0200 


0.0300 


0.1348 


0.0270 


1.4466 


0.0086 


V570_Per B 


0.5050 


0.0180 


0.0200 


0.0300 


0.1348 


0.0270 


1.3471 


0.0081 


CD_Tau A 


0.6320 


0.0160 


0.0800 


0.1500 


0.6757 


0.1351 


1.4420 


0.0160 


CD_Tau B 


0.5220 


0.0170 


0.0800 


0.1500 


0.6757 


0.1351 


1.3680 


0.0160 


AD_Boo A 


0.6400 


0.0330 


0.1000 


0.1500 


0.3387 


0.0677 


1.4136 


0.0088 


ADJBoo B 


0.2780 


0.0350 


0.1000 


0.1500 


0.3387 


0.0677 


1.2088 


0.0056 


VZ_Hya A 


0.4800 


0.0390 


-0.2000 


0.1200 


0.3387 


0.0677 


1.2713 


0.0087 


VZ_Hya B 


0.2410 


0.0420 


-0.2000 


0.1200 


0.3387 


0.0677 


1.1459 


0.0059 


V505_Pcr A 


0.4270 


0.0210 


-0.1200 


0.0300 


0.3387 


0.0677 


1.2719 


0.0072 


V505_Per B 


0.3990 


0.0210 


-0.1200 


0.0300 


0.3387 


0.0677 


1.2540 


0.0072 


UX_Men A 


0.3820 


0.0290 


0.0400 


0.1000 


0.5367 


0.1073 


1.2350 


0.0058 


UX_Men B 


0.3200 


0.0300 


0.0400 


0.1000 


0.5367 


0.1073 


1.1957 


0.0072 


ALPhe A 


0.6870 


0.0440 


-0.1400 


0.1000 


0.8507 


0.1701 


1.2336 


0.0045 


ALPhe B 


0.6720 


0.0430 


-0.1400 


0.1000 


0.8507 


0.1701 


1.1934 


0.0041 


WZ_Oph B 


0.4060 


0.0290 


-0.2700 


0.0700 


0.6757 


0.1351 


1.2268 


0.0071 


WZ_Oph A 


0.4030 


0.0290 


-0.2700 


0.0700 


0.6757 


0.1351 


1.2201 


0.0062 


V 432 Aur B 


0.4290 


0.0230 


-0.6000 


0.0500 


0.8507 


0.1701 


1.0786 


0.0053 


alf Cen A 


0.1900 


0.0080 


0.2400 


0.0400 


0.8507 


0.1701 


1.1050 


0.0070 


alf Cen B 


-0.3030 


0.0210 


0.2400 


0.0400 


0.8507 


0.1701 


0.9340 


0.0060 


NGC188JCR_V12 A 


0.3440 


0.0320 


-0.1000 


0.0900 


1.3482 


0.2696 


1.1034 


0.0074 


NGC188_KR_V12 B 


0.3050 


0.0320 


-0.1000 


0.0900 


1.3482 


0.2696 


1.0811 


0.0068 


V568JLyr A 


0.2580 


0.0320 


0.4000 


0.1000 


1.6973 


0.3395 


1.0745 


0.0077 


V568_Lyr B 


-0.5150 


0.0360 


0.4000 


0.1000 


1.6973 


0.3395 


0.8273 


0.0042 


V636 Cen A 


0.0530 


0.0250 


-0.2000 


0.0800 


0.2690 


0.0538 


1.0518 


0.0048 


V636 Cen B 


-0.4130 


0.0350 


-0.2000 


0.0800 


0.2690 


0.0538 


0.8545 


0.0030 


Sun 


0.0000 


0.0001 


0.0000 


0.0050 


1.0000 


0.0000 


1.0000 


0.0001 



7 



